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Abstract. An overview of the related topics of anomalous coarsening and glassy 
dynamics is given. In anomalous coarsening, the typical domain size of an ordered 
phase grows more slowly with time than the power law dependence that is usually 
observed, for example, in magnetic systems. We discuss how anomlaous coarsening may 
arise through domain-size dependent energy barriers in the coarsening process. We also 
review the phenomenology of glassy dynamics and discuss how simple nonequilibrium 
models may be used to reproduce certain aspects of the phenomenology. In particular, 
models involving dynamical constraints that give rise to anomalous coarsening are 
considered. Two models, the Asymmetric Constrained Ising Chain and the ABC 
model, are discussed in detail with emphasis on how the large energy barriers to 
coarsening arise through the local dynamical constraints. Finally, the relevance of 
models exhibiting anomalous coarsening to glassy systems is discussed in a wider 
context. 
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1. Introduction 

In this paper I will review some very simple dynamical models whose dynamics slow 
down with time. Along with this, some typical length scale in the system grows. 
These two features are referred to as 'glassy dynamics' and 'coarsening' respectively. 
Anomalous coarsening refers to the situation where a length scale grows more slowly 
with time than the usual power law dependence. In this introduction I will briefly review 
these topics and their interconnections. Then in the later sections I will describe in detail 
two models that exhibit anomalous coarsening and glassy dynamics. The appealing 
feature of the models to be discussed is that although they exhibit non-trivial behaviour 
they are simple enough to analyse and gain a firm understanding of. In Sections 2,3 I 
will summarise these analyses. Although the choice of models to be reviewed reflects 
personal research interests, I believe that the two models studied in Sections 2 and 3 are 
each representative of a class of systems. To bring this out I try to make connections 
to other related models in Sections 2.5 and 3.8. Finally, in Section 4, I return to the 
relation between glassy dynamics and coarsening. 

1.1. Glassy dynamics 

The phenomenolgy of glassy systems is well established — see ^ for excellent reviews. 
Experimentally, the archetypal system is a liquid that, when cooled rapidly to 
temperatures where the equilibrium state should be a crystalline solid, becomes trapped 
in metastable liquid-like configurations. Thus the higher temperature configurations are 
frozen in and one can meaningfully say that the glass is like a frozen liquicl|[]]. 

Three distinguishing features of the glassy state are the long relaxation times, 
stretched exponential decay of correlation functions (see equation 3 below) and ageing 
phenomena whereby since the system is out of thermal equilibrium it evolves 
continuously as time goes by and time-translational invariance is lacking. This 
phenomenology provides a mandate for theoretical study. 

The long relaxation times, that show a non-Arrhenius divergence as the temperature 
T is lowered are often fitted experimentally by the Vogel-Tammann-Fulcher (VTF) law 

r = roexp[-B/(T-ro)] . (1) 

The relaxation time r may characterise, for example, the time for an externally imposed 
stress to relax. Although some heuristic justifications have been offered f^, for practical 
purposes VTF is just a fit with three parameters to,B,To. Tq = reduces to an 
Arrhenius law. A system for which Tq is small, so that one has something close to 
Arrhenius behaviour, is referred to as a 'strong glass', whereas a system exhibiting large 
deviations from Arrhenius behaviour is referred to as a 'fragile glass'. Generally, Tq 
is much lower than the experimental temperatures so although there is a singularity 
in the fit, it is not physically relevant. Nevertheless it should be noted that there has 

II Strictly one may distinguish between a supercooled liquid and a glass according to the rate of the 
cooling schedule but this issue is not pertinent here. 
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been a long debate concerning whether Tq represents a true thermodynamic transition 
temperature achievable in the limit of infinitely slow cooling. 

On the other hand, other alternative functional forms for relaxation times r(T) 
have been proposed that do not exhibit singularities at any finite T. Among these, the 
exponential inverse temperature squared (EITS) form 

T ~ exp(const/T^) (2) 

(where the constant is positive) is popular. Experimentally, it is difficult to distinguish 
between VF and EITS behaviour due to obvious limitations on the longest accessible 
timescales; both can represent the experimentally observed t(T) in many materials @]. 

Stretched exponential decay of correlation functions, let us say an autocorrelation 
q{t), is expressed by the Kohlrauch- Williams- Watt law 

g(t)~exp[-(t/r)T (3) 

where the stretching exponent 6 < 1. An heuristic explanation for this law is to postulate 
a broad distribution f2(r) of relaxation modes with decay constants r yielding 

q{t) = J dr Q{r)exp{-t/T). (4) 

If one assumes ^{t) ~ exp(— ar''') then for large t a poor man's saddle point argument 
implies that the dominant modes have r = {t/a'yy^^"''^^^ which leads to @ with 
6 = 7/(1 + 7). However the question remains as to how the broad distribution of 
modes f2(r) comes about. 

1.2. Kinetically constrained models 

One idea that was proposed to generate a broad distribution of relaxation times was of 
a hierarchy of degrees of freedom f^. The different levels in the hierarchy then relax in 
serial, the degrees of freedom in one level having to wait for the degrees of freedom in 
the level above to reach some configuration before they are free to evolve. This latter 
condition is a realisation of a dynamical constraint. 

A more concrete realisation of a dynamical constraint in a system with just one 
set of degrees of freedom is the n-spin facilitated kinetic Ising models introduced by 
Fredrickson and Anderson That model comprises non-interacting Ising spins in 

a downwards pointing field. However a spin can only fiip if at least n nearest neighbour 
spins are pointing up (against the field). This gives rise to slow cooperative relaxation. 



A modification of this model is to have anisotropic dynamical constraints In 
the one-dimensional version which we will refer to as the Asymmetric Constrained Ising 
Chain (ACIC) a spin can fiip down only if its left neighbour is pointing up. As we shall 
show in Section ^, although the equilibrium distribution is the Boltzmann distribution, 
the relaxation to equilibrium is strongly affected by the asymmetric constraint. Thus 
when quenching from high temperature to low temperature the equilibrium distribution 
implies that most spins should be pointing down. However for a spin to flip down its 
neighbour has to point up. Thus there is an energy barrier for isolated up-spins to flip 



Anomalous coarsening and glassy dynamics 



4 



up that is related to the size of the domains of down-spins separating these up-spins (this 
will be quantified in section 0). As the system gets closer to equilibrium these domains 
become longer and the energy barriers increase. Thus the dynamic slows down. 

Dynamical constraints can also induce 'entropy barriers'. In this case there are no 
energy barriers to relaxation, rather one can imagine special configurations or small 
doors in the phase space that the system must pass through to allow it to relax. 
These doors are found through random exploration of the phase space. However, as 
the relaxation proceeds the doors become progressively fewer and harder to find. One 
example of such a system is the Backgammon Model [[l5l . 

Finally one should contrast the idea of energy (or free energy) barriers induced 
by dynamical constraints with energy barriers induced by disorder. It is well known 
that in a disordered system, where there is competition between different quenched 
random interactions, one can have large energy barriers in the phase space. It has been 
argued that such quenched disorder can mimic glassy systems (which are generally non- 
disordered). The theoretical machinery developed in the study of disordered systems 
then allows one to to proceed in calculations |T^. However, one still has to come up with 
arguments that relate the quenched disorder to some dynamically 'self-induced disorder' 
ir?] , |18|. Alternatively one can go one step further and simply make assumptions about 
the phase space, such as a valley structure, without specifying how such features arise as 



a result of the microscopic model. This forms the basis of the 'trap model' |T^, |2^. In 
that scenario one has a distribution of trap depths in the phase space. With increasing 
time the system will explore deeper and deeper traps and remain in them for longer and 
longer. 



1.3. Coarsening 

We now turn to the idea of coarsening and make comparisons with the glassy dynamics 
discussed so far. To visualise a coarsening system think of quenching a system from a 
high temperature phase where its order parameter is zero to a low temperature where 
the order parameter can take some number of non-zero values (each different value 
corresponding to a different ordered phase). Domains of the ordered phase(s) emerge 
and grow in time and it is this phenomenon that is referred to as coarsening. At late 
times the system enters a scaling regime, that is a regime characterised by a single length 
scale (the typical size of domains) i{t) that grows as i{t) ~ t". In this late time scaling 
regime the distribution of domains, once scaled by £{t), is statistically invariant. Thus 
the typical domain size indicates the age of the coarsening system. The value of the 
exponent n depends on the symmetry of the order parameter and conservation laws of 
the system. For a review see ||22|| . 

As a concrete example consider the zero temperature Ising model with Glauber 
spin-flip dynamics (non-conserved order parameter). Domains of up-spins and down 
spins are separated by domain walls that perform random walks: a step of the walk 
corresponds to the event that one of the spins adjacent to the domain wall flips. When 
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two domain walls meet they annihilate and a domain is eliminated. Straight away one 
can argue roughly that, since the domain wall motion is diffusive, the time typically 
required to eliminate domains of size / is T{1) ~ Thus the typical domain size after 
time t is £ ~ t^/^. 

More precisely one can write the growth law as a differential equation by noting 
that the rate of change of the typical domain size should be proportional to the inverse 
of the mean time to eliminate a domain multiplied by the size of the domain being 
eliminated. 

and one recovers i ~ t^^"^. This growth exponent actually holds for a non-conserved 
scalar order parameter in two dimensions and above 



Now consider generalising (|^) to the case where some energy barrier AE (or more 
strictly free energy barrier) is involved in the elimination of the domains and the system 
is at low but finite temperature. 

Actually this precludes one dimensional systems that only order up to a finite length 
scale at finite temperature (see Section but at very low T one can consider the 
ordering process up to that finite length. Some possible scenarios resulting from (|^) 
have been categorised in |2T|. If the barriers IS.E are independent of £ one recovers 
i ~ t^/^ growth. If the barriers are proportional to one obtains i ~ [Int]^'''", in 
particular m = 1 yields logarithmic growth. In Section ^ it will be shown that for the 
ACIC discussed in Section |1.2| the energy barriers are logarithmic in I thus yielding 
from (^ a growth law where the growth exponent is proportional to the temperature as 
T — s> 0. In section |] models will be discussed that have energy barriers that are linear in 
I thus yielding domain growth that is logarithmic in time. For a particular model, the 
ABC model [23| it will be shown how these linear energy barriers arise. We refer to such 
cases where something different from power law growth with temperature independent 
exponent is exhibited as anomalous coarsening. 

A common approach in the study of coarsening is to make an approximation 
of a mean-field nature. That is, one focusses on the probability distribution of the 
domain sizes and ignores spatial correlations between domains. Such an approach is 
variously referred to as an Interparticle Distribution Function |2|], Independent Intervals 
Approximation p5[ or a 'bag model' p6|. To visualise this one thinks of placing all the 
domains in a bag i.e. forgetting how domains are arranged with respect to one another. 
Then the dynamics become the updating of the domains in the bag. At each update a 
domain is selected from the bag along with temporary neighbours. An update rule is 
implemented that depends on the shape/ size of the domain and its neighbours. After the 
update the domains are all replaced in the bag. For example, in models where domain 
walls diffuse and annihilate or coalesce, such as the kinetic Ising model discussed above, 
the update rule is to lengthen and shorten domains according to how the domain wall 



Anomalous coarsening and glassy dynamics 



6 



moves. In another class of models 0, domains are eliminated at each update 
with probabilities depending on their size; then the length of an eliminated domain is 
distributed amongst the neighbours. 



1.4- Nonequilibrium steady states 

So far, although we have discussed glassy systems out of thermal equilibrium, the 
dynamics have implicitly been assumed to obey detailed balance in the equilibrium 
state. Detailed balance means that in the steady state there is no net flow of probability 
between any two configurations. However one can consider a much larger class of 
nonequilibrium systems that are defined solely by their dynamics, without reference 
to any energy function. Although these systems may relax to some steady state it 
need not be a steady state described by Gibbs-Boltzmann statistical mechanics. In 
general there will be a net flow of probability between pairs of configuration, leading to 
probability current loops in the configuration space. 

Examples of nonequilibrium steady states are given by driven systems with open 
boundaries where a mass current is driven through the system. Thus the system is 
driven by its environment rather being in thermal equilibrium with its environment. In 



such a driven steady state generic long range correlations may be exhibited |2^. This 
is in contrast to an equilibrium state which only exhibits long-range correlations at 
non-generic points i.e. phase transitions. Since there is no energy function and Gibbs- 
Boltzmann statistical mechanics does not apply, there is no general formulation within 
which to solve for such a driven steady state. However it has often been suggested 
that the generic long-range correlations may result from some effective long-range 
Hamiltonian that could describe the dynamics. 

Of particular interest have been one- dimensional nonequilibrium systems. For 
models respecting detailed balance it is well known that no phase transition or ordering 
process that continues indefinitely can occur. However for nonequilibrium systems this 



is not the case 1 29]. Thus nonequilibrium systems afford new possibilities for coarsening 



processes even in one-dimensional systemsl^O]]. Moreover there are a number of exactly 



solvable one dimensional nonequilibrium systems |31[]. In section ^ the steady state of the 
ABC model will be solved for some special cases and it will be shown how strong phase 
separation along with anomalous coarsening can occur in one dimensional systems. 



2. Asymmetric Constrained Ising Chain 

In this section I discuss the model introduced by Jackie and Eisinger|P, |lOl. As discussed 
in Section |1.2| it was originally introduced as a model of cooperative, glassy relaxation. 
In particular the directed nature of the dynamical constraint implies a hierarchy in the 
spin relaxation. As I shall now describe the directed nature of the constraint also makes 
the model amenable to analysis. 
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2.1. Model definition 

The model comprises L Ising spins Sj = 0, 1 on a one- dimensional lattice with periodic 
boundary conditions (site i = L + 1 is identified with site i = 1). The dynamics are 
defined by the following spin-flip rates 

11^10 with rate 1 

(7) 



1 — i> 1 1 with rate e 



where 



e = exp(-l/T) . (8) 

Thus a spin can only flip if its left neighbour is pointing up (note that in ||^ the mirror 
image of the above deflnition was used, so that that right neighbour had to point up 
for a spin to be able to flip). By a rate, say x, we mean that in a small time dt the 
event happens with probability xdt. It is easy to check that the dynamics obey detailed 
balance with respect to an energy function E = Yli=i ^- ^- equilibrium distribution 
corresponds to free spins in a downwards pointing fleld: 



Pcq{{si}) = ^exp 



(9) 



T 

It follows that the equilibrium concentration of up-spins, c = (si), is given by 
e 



1 + e)^ ■ 

Dins, c = (si), is given by 

(10) 

1 + e ^ ^ 

Although (^) would hold for a number of different dynamics obeying detailed balance, 
other properties such as two time correlation functions may be sensitive to the particular 
choice of dynamics. 

It is clear that the asymmetric dynamical constraint implies that information 
propagates to the right only. Thus in a thermodynamic limit where information 
cannot propagate all the way around the ring back to the starting point, we must have 
(si+fc(0)sj(t))eq — = for /c > 0. However, the fact that detailed balance holds implies 
that in the steady state we must have reversibility. To see this note that when one 
has detailed balance there is no net flow of probability between any two conflgurations. 
Since there is no flow of probability there is nothing to distinguish the forwards direction 
in time from backwards direction. Therefore running the systems backwards in time will 
not change any two time correlation functions and 

{Si+k{0)Si{t))eq = {Si+k{t)Si{0))cci . (11) 

Since we have argued that for A; > the left hand side of (|ll]) is equal to we deduce 



that the connected correlation function must be site diagonal: 

{s,{0)s,{t)),^-c^ ^6ij . (12) 

This result is particular to the fully asymmetrically constrained model. 

We will be interested mainly in the behaviour after a quench from equilibrium 
at some high initial temperature T 1 to low temperature T -C 1 (e — >^ 0). At 



Anomalous coarsening and glassy dynamics 



8 



T>1 0110101110011 
T< 1 0100101000010 
c^e/(l + e) 0100000000010 



Figure 1. Schematic representation of the evolution of the system following a deep 
quench. Before the quench T ^ 1 and c, the concentration of up-spins, ~ 1/2. 
After the quench all the 'mobile' up-spins (i.e. those adjacent to other up-spins) are 
eliminated first. A slow coarsening process ensues that reduces the density of up-spins 
c to its equilibrium value. 



low temperatures the equilibrium concentration of up-spins (0) is small. Thus the 
quench is followed by a process of elimination of up-spins. However to eliminate an up- 
spin one first has to generate an adjacent up-spin. This implies energy barriers in the 
system's evolution. In figure the sequence of events after such a quench is illustrated 
schematically. 

The basic objects that we use for the description of the system are domains. As 
shown by the vertical lines in 

...1|0001|1|1|01|001|1|1|01|0..., 

a domain consists of an up-spin and all the down-spins that separate it from the nearest 
up-spin to the left. The length d of a domain then gives the distance between the up- 
spin at its right edge and the nearest up-spin to the left. Note that adjacent up-spins 
are counted as separate domains of length d = 1. In equilibrium, the distribution of 
domain lengths and its average are 

Peq(rf) = e/(l + rfeq = 1 + 1/e • (13) 

Now consider what happens after a deep quench to T ^ 1, e ^ 1. The equilibrium 
concentration of up-spins at the final temperature T is c = 1/deq = e + 0{e^)] hence the 
equilibrium probability of finding an up-spin within a chain segment of finite length d 
is 0{de) and tends to zero for e ^ 0. In this limit (e ^ at fixed d), the flipping down 
of up-spins therefore becomes irreversible to leading order. In terms of domains, this 
means that the coarsening dynamics of the system is one of coalescence of domains: an 
up-spin that flips down merges two neighbouring domains into one large domain. 

Such coarsening processes have been studied in a variety of contexts. In particular 
irreversible coarsening processes in which the rate of elimination depends solely on the 
domain size have a very convenient property: during such a process, no correlations 
between the lengths of neighbouring domains can build up if there are none in the 
initial state [|^. For the present model the equilibrated initial state consists of domains 
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d = 2 101^ 111 110^ 100 h{2) = 2 
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d = A 10001 ^ ^11101^10101^ 

lOlii^ ^10100^ ^ 10000 /i(4) = 3 



Figure 2. Paths through spin configurations in the ehmination of a domain of size d 
that traverse the minimum energy barrier. The height of the barrier is h{d) and the 
highest energy configuration(s) along the path has its excess excited spins underhned 



independently distributed according to (|T3D and is indeed uncorrelated. We take 
advantage of this property in Section |2.3| to obtain an exact solution of the coarsening 
dynamics. We first discuss in detail how energy barriers arise in the dynamics. 

2.2. Energy barriers 

First we estimate the typical rate T{d) at which domains of length d disappear by 
coalescing with their right neighbours. Because domain coalescence corresponds to the 
flipping down of up-spins, r{d) can also be defined as follows. Consider an open spin 
chain of length d, with a 'clamped' up-spin (sq = 1) added on the left. Starting from the 
state (so, si, . . . ,Sd) = 10 . . . 01, T~^{d) is the typical time needed to reach the empty 
state 10 ... 00 where spin Sd has 'relaxed' i. e. has flipped down. Any instance of this 
relaxation process can be thought of as a path connecting the initial and final states. Let 
us call the maximum number of 'excited' spins (up-spins except sq) encountered along 
a path its height h. One might think that the relaxation of spin Sd needs to proceed 
via the state 11. . . 1, giving a path of height d. In fact, the minimal path height h{d) is 
much lower and given by 

h{d) =n+l for 2"~^ < d < 2" (14) 

where n = 0,1, . . .. 

To get a feeling for the result ([T^) consider in figure § some small domain sizes. 
The figure illustrates that to generate an up-spin adjacent to the right boundary spin 
of the domain one can proceed via a sequence of stepping-stones e.g. for d = 4 one first 
generates an isolated up-spin in the middle of the domain then uses this stepping-stone 
to generate the subsequent excited spins in a similar manner to the relaxation of a = 2 
domain. 

The result (|l^ is easily demonstrated for = 2" [|13|. To relax the 2^-th spin 



52", one can first flip up S2n-i and use it as a stepping-stone for relaxing 52". The 
corresponding path is (with 82^-1 and S2" underlined) 
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and reaches height /i(2") = /i(2"~^) + 1; the +1 arises because the stepping-stone stays 
up while the spin 2^~^ to its right is relaxed. Continuing recursively, one arrives at 
/i(2") = h{l) + n; but h{l) = 1 because the only path for the relaxation of si is 
11 — s> 10. Thus we obtain equation (|T4|) for d = 2"^; a proof for general integer d is 
given in [Q. 

From (|T^ it is evident that the energy barrier AE for the elimination of a domain 
of size d is AE ~ Ind/ In 2. Thus the rate at which such domains are eliminated is 

From the discussion of Section |L3| and equation ^ we deduce that the typical domain 
size grows and the typical energy (number of up-spins) decreases as 

E,,p~r^'"^ (16) 

Also since d^^ — e = e^^'^ the equilibration time is 

teq~exp [1/T2ln2] . (17) 



2.3. Hierarchical coarsening 

From the scaling of T{d), the coarsening dynamics in the limit e — > naturally divides 
into stages distinguished by n = h{d) — 1 = 0, 1, . . .. During stage n, the domains 
with lengths 2"~^ < (i < 2" disappear; we call these the 'active' domains. This process 
takes place on a timescale of 0(T~^{d)) = C(e~"); because the timescales for different 
stages differ by factors of 1/e, we can treat them separately in the limit e ^ 0. Thus 
during stage n active domains are eliminated and the distribution of inactive domains 
{d > 2"^) changes because elimination of an active domain implies coalescence with a 
neighbouring domain and results in the creation of a new inactive domain. 

As discussed above for the irreversible system there are no correlations between 
neighbouring domains. Therefore we can work directly with the probability of domain 
sizes P{d, t) i.e. the Independent Intervals Approximation sketched in Sec is actually 
exact. 

First let us write down the master equation during stage n where we assume that 
all domains with d < 2"~^ have been eliminated. 

— P(rf, t) = -T{d)P{d, t)+ ^{d')P{d\ t)P{d - d', t) . (18) 

d'=2"-'-+l 



The first term in the right hand side of ([TsD represents domains of size d being eliminated; 
the second term represents the domains of size d being create through a domain of size 
d' coalescing with a domain of size d — d'. Now we introduce rescaled time r = te"; 
during stage n of the dynamics and in the limit e— >0, it can take on any positive value 
r > 0. Defining T{d) = T{d)/e^ and taking the limit e — > the master equation reduces 
to 

for 2''>d> 2"-^ ^) = -f (rf)P(rf, r) (19) 
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on 

for rf>2" p(d,T)= V r(d')P(d',t)P(d-d',t) . (20) 

To proceed we define the generating function 

oo 

G{z,t)= Y1 P{d,r)z' (21) 
and its analogue for the active domains, 

2" 

H{z,r)= nd,r)z'. (22) 

d=2"-i+l 

Then multiplying (^,^) by z'^ and summing appropriately yields 

^H{z,t)= - V{d)P{d,T)z'' (23) 



2"-i+l 

oo 2" 



G(z,r)= i7(^,r)+ ^ ^ f(ciOP(ci', r)P(ci - d', r)/ 



dr ' dr 

rf=2"+l d'=2"-l+l 



(^^ff(rf,r))[l-G(2,r)] (24) 



(where the last equality follows by reordering the sums over d and d'). Equation |2^ may 
be integrated and one obtains 

[Zct^l] = - [Hi^, r) - H{z, 0)] . (25) 

Now at the end of stage n, all domains that were active during that stage have 
disappeared, and so H{z, oo) = 0. Thus 

G{z, oo) - 1 = [G{z, 0) - 1] exp[H{z, 0)] . (26) 

Recall that we are considering stage n of the dynamics. The initial condition for 
stage + 1 of the dynamics will be given by the distribution P{d, t) at the end of stage 
n. Thus defining G„ = G{z, 0) for stage n, with a similar definition for the active 
generating function we can relate the different stages of the dynamics through 

Gn+i{z) - 1 = [Gn{z) - 1] exp[Hr,{z)] . (27) 

This exact result relates, through their generating functions, Pn{d) and Pn+i{d) which 
are defined as the domain length distributions at the end of stages n — 1 and n of the 
dynamics respectively. Iterating it from a given initial distribution Po{d) gives Pn{d) 
for all n = 1,2, . . .. Figure ^ shows numerical results for the case where Po{d) is the 
equilibrium distribution ([T3|) corresponding to an initial temperature of T = oo. It 
is clear that a scaling limit emerges for large n. By this it is meant that rescaled 
distributions 

P„(x) = 2"-ip„(d) where x = , (28) 




Figure 3. Domain length distributions Pn{d) at the end of stage n — 1 of the low 
T coarsening dynamics, for initial temperature T = oo. Open symbols and lines: 
Theoretical results, calculated from (p7|), for n = (Q; initial condition), 1 (Boxes), 
2 (diamonds), 3 (triangles). Full symbols: Simulation results for a chain of length 
L = and e = lO^"* {n = 1,2) and e = lO'^ (n = 3). Inset: Scaled predictions 
2"~^Pn{d = 2"^^x) vs. X for n = 1, . . . , 8. Bold line: Predicted scaling function (|33|). 
(Figure taken from p4|). 

converge to a limiting distribution P{x) for the scaled domain size x. This is just a 
statement of the invariance of the coarsening processes in each stage once the domain 
sizes are rescaled by the characteristic size domain size 2"~^. 

The change to a continuous variable x for the domain lengths simply results in 
generating functions G{z,t), H{z,t) being replaced by Laplace transforms. Invariance 
under (^) then gives the equation 

gi2s)-l = [g{s)-l]exp[h{s)] (29) 

where 

dxP{x)e~''' h{s,T) = j dxP(x)e-'^. (30) 

We found a solution to this equation by noting that the numerics strongly suggest 
P{x) = 1/x for 1 < X < 2. Using this as an ansatz implies 



h{x) = Ei(s) - Ei(2s) where Ei(s) = 




which when inserted into (^) yields 

[1 — g{s)] exp(Ei(s)) = constant (31) 
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The requirement that g{s) — > for large s fixes the constant as unity which yields 

g{s) = l-exp{-Ei{s)). (32) 

Expanding the exponential as a series allows the Laplace transform to be inverted term 
by term and one obtains 

Pix) = — / e^^g{s) 

1 /•■y+ioo °° 17';'^ 

2m Jj-ioc m\ 



e(x-i)--e(x-2) ^'^^^-^^ + (33) 

X X 



where G(x) is the Heaviside step function. This series ( p3|) has singularities in the /c-th 
derivative at the integer values a; = fc + l, fc + 2, ... 

It is interesting to note that P{x) given by (|33D is identical to the scaling function 
obtained for a simple 'paste-all' model of coarsening wherein the smallest domain on a 



one dimensional lattice is eliminated by pasting it onto one of its neighbours [p6[]. 

The calculated P{x) agrees well with the results obtained by direct iteration of (p7|) 
(figure The average domain length in the scaling limit is given by = 2"~^x; from 
the results for P{x) we find x = exp(7) = 1.78 . . ., where 7 is Euler's constant. 

2.4- Stretched exponential relaxation 

The result (0) for the EITS equilibration time teq = exp(l/T^ln2) is based on the 
extrapolation of the finite-d coarsening behaviour, ~ t^^"^^, into the equilibrium region 
d = deq = C(l/e), where it is no longer strictly valid. We now show, however, that 
the same timescale is obtained from the initial decay of the spin-spin autocorrelation 
function defined by 

= M5M)k (34) 

(Si(0))cq 

at equilibrium at low temperature T. 

In equation |3^ R{t) is the probability that an up-spin at t = is also up at a 
later time t. As t increases, R{t) decays from 1 to the equilibrium concentration of 
up-spins, c = e/(l + e). To find the initial decay of R(t), consider the early stages of the 
dynamics {i.e. t = 0{e^'^) with u finite). For u n+, all domains of length c? < 2" will 
have disappeared because t ^ T^^{d). Therefore only up-spins that bounded longer 
domains at t = will have an 0{1) probability of still being up. From the equilibrium 
distribution (|I^), one sees that they constitute a fraction (1 + e)~^" of the up-spins at 
t = 0, hence R{iy = n -|- 0) ~ 1 — 2"e -|- O(e^). Neglecting corrections of C(e^), the 
quantity — lni?(z/) thus lies between and (for u > 0). 
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Reverting to unsealed time t, we have 

l/2<-[lni?(t)]/(t/teqri°2<i (35) 

for short times (t/teq)^^°^ ^ 1 which imphes 



R ~ exp — 



a(tAeq)''"^'j (36) 
i.e. stretched exponential relaxation with a stretching exponent that depends linearly on 
T. However this argument only holds strictly for (t/teq)"^'"^^ <^ 1- For longer times the 
stretched exponential behaviour no longer holds and one requires a more sophisticated 
analysis P2[] . 

Note that the timescale teq that enters here is the same as the equilibration time 
teq = exp(l/T^ ln2) found above. Thus to leading orderwe can identify the equilibration 
time for coarsening after a quench, with the equilibrium relaxation time; both have an 
EITS-divergence at low T. However the corrections to the EITS-divergence for the two 
timescales ( e.g. factors of the form exp(a/T)) need not be equal. 

2.5. Other related models 

Recently it has been shown that the coarsening theory described above for ACIC 



is also relevant to a three-spin interaction Ising Model on a triangular lattice |33 
There a dual model entails the elimination of defects (corresponding to up-spins in 
the ACIC), subject to dynamical constraints. An EITS relaxation time is obtained and 
the Independent Intervals Approximation described here (exact for ACIC) can be used 
as a good approximation to the coarsening dynamics (see Garrahan this volume). 

One can also interpolate between the asymmetric constrained model and the 



symmetric constrained model by introducing a parameter b into the dynamics |34 
The spin-flip rates are 

1-6 

11 ^ 10 

{l-b)e (37) 

11 ^ 1 

be 

The ACIC is recovered when 6 = and the symmetric model is recovered when 
6 = 1/2. In the symmetric model an isolated up-spin can effectively diffuse by creating a 
neighbouring up-spin then flipping down the original up-spin. The domain coalescence 
happens through this diffusion process. Thus the energy barrier for elimination of 
domains does not depend on domain size. This implies an Arrhenius relaxation law 
and 'strong' glass behaviour. The crossover to the 'fragile' glass behaviour seen for the 
ACIC {i.e. EITS relaxation time) as one varies the parameter 6 has been studied (see 
Buhot this volume). 

Finally let us mention a constrained Ising spin chain where the field is induced 
dynamically The allowed spin flips have rates defined as follows 
1/2 

011^:^001 

1/2 
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110^ 100 (38) 
1/2 

010^ 00 

Thus a down-spin inside an up-spin domain cannot flip up but an up-spin inside a down- 
spin domain can flip down. The domains of up-spins grow as normal as t^/^ but the 
down-spins domains grow slightly more quickly as t^^^ In t. Eventually this results in a 
slow decay of the magnetisation (number of up-spins as) c ~ 1/lnt. 



3. The ABC Model 



3.1. Coarsening in one- dimensional systems 

First let us review why indefinite coarsening does not occur in equilibrium systems in 



one dimension. Perhaps the best known argument is that of Landau and Lifshitz ||36 
(Other arguments are summarised in [Q.) For simplicity, consider a one- dimensional 
lattice of L sites with two possible states, say A and B, for each site variable. Let us 
assume the ordered phases, where all sites take state A or all sites take state B, have the 
lowest energy, and assume a domain wall (a bond on the lattice which divides a region 
of A phase from that of B) costs a finite amount of energy e. Then n domain walls will 
cost energy ne but the entropic contribution to the free energy due to the number or 
ways of placing n walls on L sites ~ nT [ln(?7,/L) — 1] for 1 -C tt, <C L. Thus for any 
finite temperature a balance between energy and entropy ensures that the number of 
domain walls grows until it scales as L, that is, until the typical ordered domain size is 
finite. 

Note that this argument relies on a finite energy cost for domain walls, and short- 
range interactions so that one may ignore the interaction energy of domain walls. If the 
domain walls can feel each other through some long-range mechanism then coarsening 
can ensue |37]. Also, of course, we require non-zero temperature so that entropy comes 
into play. In contrast, at zero temperature the one-dimensional kinetic Ising model 
discussed in the introduction does coarsen. 

In the following we will discuss a one-dimensional model where although the 
dynamics are local the systems coarsens. In a special case one can understand this 
through the existence of an effective long-range energy function. 



3.2. Model definition 

Here we define a model, to be referred to as the ABC model, that exhibits phase 
separation in one dimension. Consider a one-dimensional periodic lattice of length 
N where each site is occupied by one of the three types of particles. A, B, or C. The 
model evolves under a random sequential update procedure which is defined as follows: 
at each time step a pair of neighbouring sites is chosen randomly and the particles at 
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these sites are exchanged according to the following rates 

q 

AB ^ BA 
1 

BC ^ CB (39) 
1 

q 

CA ^ AC. 
1 

The particles thus diffuse asymmetrically around the ring. The dynamics conserve the 
number of particles, Na, Nb and Nc of the three species. 

The g = 1 case is special. Here the diffusion is symmetric and every local exchange 
of particles takes place with the same rate as the reverse move. The system trivially 
obeys detailed balance reaching a steady state in which all microscopic configurations 
(compatible with the number of particles A'^, Nb and Nc) are equally probable. This 
state is disordered and homogeneous; no phase separation takes place. 

a) C ABAC BC AABBC BC A 
h) AAABBCC AABBBCCC 
c) CC AAAAABBBBBCCC 



Figure 4. Schematic representation of the evolution of the system starting from a 
random initial condition a). Initially all unstable domain walls {B A, C B or AC) are 
removed and one arrives at a metastable state b) . A slow coarsening process ensues in 
which the smallest domains are eliminated until one arrives at a fully phase separated 
state c). Compare with figure |l|. 

Now consider the case g < 1 (the case q> 1 can easily be understood by symmetry). 
As a result of the bias in the exchange rates an A particle moves preferentially to the left 
inside a B domain and to the right inside a C domain. Similarly the motion of B and C 
particles in foreign domains is biased. Consider the dynamics starting from a random 
initial configuration figure ^ a). The configuration is composed of a random sequence of 
domains of A, B, and C particles. Due to the bias a local configuration in which an A 
domain is placed to the right of a -B domain is unstable and the two domains exchange 
places on a relatively short time scale which is linear in the domain size. Similarly, AC 
and CB domains are unstable too. On the other hand AB, BC and CA configurations 
are stable and long-lived. Thus after a relatively short time the system reaches a state 
of the type illustrated in figure ^ b) in which A, B and C domains are located to the 
right of C, A and B domains, respectively. The evolution of this state takes place via a 
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slow diffusion process in which, for example A particles have to diffuse against a drift 
over an adjacent B domain. The time scale for an A particle to cross is oc q~\ where 
I is the size of the B domain. We use the discussion of Section |L3| and equation ^ to 
deduce that the system coarsens with an average domain size that increases with time 
as lnt/|lng|. Eventually the system phase separates into three domains of the three 
species of the form A . . . AB . . . BC . . .C. 

In a finite system the phase-separated state may further evolve and become 
disordered due to fluctuations. However, the time scale for this to happen grows 
exponentially with the system size. For example it would take a time of order of 
q-mm{NB,Nc} f^j, ^]-^g ^ domaiu in the totally phase separated state to break up into 
smaller domains. Hence in the thermodynamic limit, this time scale diverges and the 
phase separated state remains stable provided the density of each species is non-zero. 
Note that there are always small fluctuations about a totally phase separated state. 
However, these fluctuations affect the densities only near the domain boundaries. They 
result in a finite width for the domain walls (the density profile is not a step function 
but is smeared out like a Fermi function). The fact that any phase-separated state is 
stable for a time exponentially long in the system size amounts to a breaking of the 
translational symmetry i. e. there are N equivalent ground states but the system has to 
spontaneously choose one of them. 

Since the exchange rates are asymmetric, the system generically supports a particle 
current in the steady state which implies that detailed balance does not hold. To see 
this, consider the A domain in the phase separated state. An A particle near the 
. . . AB . . . boundary can traverse the entire B domain to the right with an effective rate 
proportional to . Once it crosses the B domain it will move through the C domain 
with speed 1—q. Similarly an A particle near the . . . CA . . . boundary can traverse the 
entire C domain to the left with an effective rate proportional to q^^ . Once the domain 
is crossed it moves through the B domain with speed 1—q. Hence the net A particle 
current is of the order of — . Since this current is exponentially small in system 
size, it vanishes in the thermodynamic limit. For the case of A^^ = Nb = Nc, this 
argument suggests that the current is strictly zero for any A^. 

The arguments presented above suggesting phase separation for g < 1 may be easily 
extended to g > 1. In this case, however, the phase separated state is BAC rather than 
ABC. This may be seen by noting that the dynamical rules are invariant under the 
transformation g ^ 1/g together with A ^ B. 



3.3. Special Case Na = Nb = Nc 

The general argument presented in the previous subsection suggests that for the special 
case A"^ = A'^ = Nc, in the steady state {i.e. after the coarsening process) there are no 
currents for any system size. We demonstrate this explicitly by showing that the local 
dynamics of the model satisfies detailed balance with respect to a long-range asymmetric 
energy function Ti. 
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We define the occupation variables as follows: 

^ I 1 if site i is occupied by an A particle 
* I otherwise. 

The variables Bi and Ci are defined similarly. Clearly the relation Ai + Bi + Ci = 1 
is satisfied. It turns out that for the case Na = Nb = Nc = N/3 the steady-state 
distribution WN{{Xi}) corresponding to the dynamics ( PU]) may be written in terms of 
an energy function 7i: 

Wn{{X,}) = . (41) 

N~l N 

n{{x,}) = Y.Y1 + ^^^i + - • (42) 

1=1 j=i+l 

Here is the partition sum given by ^ q^^d^i})^ where the sum is over all configurations 
in which A^^ = Nb = Nc- Note that although the system is periodic and site 1 is not 
in any way special ( |^ ) appears to single out site 1. Thus it is not clear that ( ^2]) is 
translationally invariant under relabelling of the spins. 

In order to turn equation (^Tj) into a usual Boltzmann form one could define g as a 
temperature variable with 

A;T = -l/lng. (43) 

Thus, g — > 1 is the infinite-temperature limit, corresponding to the disorder state where 
each configuration is equally likely. The proof of equations (ffT],H2|) is straightforward. 



This is done by considering a nearest-neighbour particle exchange and verifying that 
detailed balance is satisfied with respect to (P^. Then we just have to check that the 



energy function is translationally invariant. We defer the proof to Section |3.7| where we 
consider a more general m-species model. 

Before proceeding further to evaluate the partition sum associated with the energy 
function (^) let us make a few observations. The ground state of the energy function 
is given by the fully separated state A . . . AB . . . BC . . .C and its translationally related 
states. It easy to check that the ground state with the A domain beginning at site 1 
has zero energy since the contribution to the sum in (|^) (coming from the A^ Cj term) 
is equal to (iV/3)^. A simple way of evaluating the energy of an arbitrary configuration 
is obtained by noting that nearest-neighbour exchanges AB —>■ BA, BC —yCB and 
CA —>■ AC cost one unit of energy each while the reverse exchanges result in an energy 
gain of one unit. The energy of an arbitrary configuration may thus be evaluated 
by starting with the ground state and performing nearest neighbour exchanges until 
the configuration is reached, keeping track of the energy changes at each step of the 
way. The highest energy is and it corresponds to the totally phase separated 

configuration A . . . AC . . . CB . . . B and its translations. Note that the majority of 
configurations have energy proportional to A^^. In Section 1X1 it will be shown that this 



implies that only the ground states and low energy excitations about them contribute 
to the equilibrium state. 
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To write 7i in a manifestly translationally invariant form we define Tij^dXi}) as 
the energy function in which site io is the origin. Namely, 

N+io-2 N+io-1 

nMX,})= J2 E [C^B,+AC,+B,A,]-{N/3)\ (44) 



where the summation over i and j is modulo N. Summing (|^) over all io and dividing 
by A^, one obtains, 

N N-1 , 

n{{X,}) = 5^ 5^(1 - + + B,A,+k) - (iV/3)2 , (45) 

1=1 k=l 

where in the summation the value of the site index {i + k) is modulo A^. In the energy 
function ( |^ ) the interaction is linear in the distance between the particles, and thus is 
long ranged. The distance is measured in a preferred direction from site i to site i + k. 
Thus the interaction is asymmetric. 



3.4- Ground States and Metastable States 



A minimum of the energy (42) is realised by a configuration with no unstable domain 
walls {B A, C B OT A C) so that any exchange of nearest neighbour particles results in an 
increase in the energy. As well as the A^ ground states there are many metastable states. 
Any metastable state is composed of a sequence of domains separated by AB^ BC and 
CA domain walls i.e. A,B and C domains follow C,A and B domains, respectively 
(see figure ^ b). Therefore each metastable state has an equal number of domains 
of each type. We shall refer to any metastable state with s domains of each type, 
with s = l,...,N/3, as an s-state; the total number of domains in an s-state is 3s. 
The s = 1 case corresponds to the ground state while s = N/3 corresponds to the 
ABC ABC . . . ABC state, composed of a total of A^ domains each of length 1. Note 
that in general the domains of an s-state need not be of equal length. 

In the coarsening process it is these metastable states that control the dynamics. 
We discuss here some properties of the states such as their number and energy. 

To obtain a bound for A/'(s), the number of metstable states with s domains of 
each species, note that the number of ways of dividing A^/3 A particles into s domains 
is (^g^^^) . The number of ways of combining s divisions of each of the three types 

r 1 3 

of particles is clearly • There are at most A^ ways of placing this string of 

domains on a lattice to obtain a metastable state (the number of ways need not be equal 
to A^ since the string may possess some translational symmetry). One therefore has 



N/3 
s — 



1 



1 



<7V(s) < A^ 



A^/3 
s — 



(46) 



Thus, the total number of metastable states is exponential in A^. 

We now consider the energy of the metastable states. It is easy to convince oneself 
that among all s-states, none has energy lower than the configuration illustrated in 
figure |. 
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Figure 5. Metastable state that has the lowest energy for a given number s of 
domains of each species. The 3(s — 1) rightmost domains are of size 1 and the three 
leftmost domains are of size {N/i — s + 1) each. 



The energy of this state, Es satisfies the following recursion relation 

= Es + N/?,-s (47) 

with El = 0. To see this from figure ^, note that the s + 1-state may be created from 
the s-state by first moving a B particle from the leftmost B domain across (A^/3 — s) 
C particles to the right, costing (A^/3 — s) units of energy. Then move an A particle 
from the leftmost A domain to the right across the adjacent B and C domains; no net 
energy change results. Thus the total energy cost of the moves is (A^/3 — s), yielding 
(^). The recursion relation (^Tj), together with i?i = 0, is readily solved to give 

B,.„-l)f-£(i^. (48) 



The energy of all metastable s-states is larger or equal to Eg as given by equation 
Note that Es increases with s. Furthermore, for finite s the energy is linear in whereas 
for s (X N the energy becomes quadratic in A^. 

Now let us consider the contribution of the metastable states to the partition sum. 
Multiplying the upper bound on the number of s-state ( ^Bj ) with the lower bound on 
the energy (^5]) one obtains an upper bound on the contribution to the partition sum 
that vanishes in the thermodynamic limit 



N 



q{s'm/-i~sis-i)/2 ^ Q fors>l, g<landAr^oo. (49) 

s-1 J\ 

However, the contribution from the ground states s = 1 is A^. Thus even though 
metastable states dominate the dynamics, they do not contribute to the partition sum 
since in (ESI) the energy grows more strongly than the entropic contribution. 



3.5. Partition Sum 

We now analyse in more detail the behaviour of the partition sum. In principle one 
wants to compute 

Zj, = J2qn{x,}) (50) 

where the sum is over all configurations in which A^^ = A^^ = Nq- First note that any 
configuration that contains unstable domain walls {i.e. is not a metastable state) can 
be associated with a metastable state by a path of decreasing energy comprising nearest 
neighbour exchange eliminating of the unstable walls. Conversely the sum over all 
configurations may be implemented by summing over all ground states and metastable 
states and the excitations about those states. It is not hard to believe that as the 
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Figure 6. Excitation about a ground state at the AB domain wall of energy E = ?>: 
one A particle has penetrated two exchanges into the B domain and a second A particle 
one exchange. This corresponds to the partition 3 = 2+1. The domain wall regions 
are underlined. 



metastable states make vanishing contributions to so do excitations about them. 
This is proven rigorously in In the following we just consider the excitations about 
the ground states. Consider figure || where a low energy excitation about a ground state 
is illustrated. The excitation is localised near the A B domain wall and comprises one or 
more A particles penetrating into the B domain (equivalently B particles penetrating 
into the A domain). The energy cost is given by the sum of the distances each A particle 
has penetrated into the B domain. Thus the total number of excitations of energy m 
at the boundary is the number of ways of dividing m into an ordered set of integers 
corresponding to the distance the first A has moved, the distance the second A has 
moved and so on. This is equal to P{m) the number of partitions of the integer m. 

Several results concerning partitions are known [p^] . First the generating function 
of P{m) is given by 

oo ^ 

J2 Q'^PM = p- (52) 

m=0 K^joo 

where (g)oo = lim (1 - g)(l - g^) . . . (i _ g") . (53) 

n— >oo 

Although a simple explicit formula for P{m) does not exist the asymptotic behaviour 
is given by 

P(m) ^ — ^ exp (7r(2/3)^/2 ^1/2^ _ (54) 
4m V 3 

Note that the increase is a stretched exponential in m i.e. slower than exponential. 

In the thermodynamic limit one can use (|52D directly to calculate the sum of 
excitations around a domain wall i.e. the sum over m of q"^, the weight of an excitation 
of energy m, multiplied by P{m), the number of such excitations. For a finite system 
there should be some upper limit on m — for example, an A particle moved across the B 
domain will eventually reach the C domain — but this upper limit can be safely taken to 
infinity for large |^|. By the same token the three domain walls have no significant 
interaction. Then one has that in the large N limit and for all g < 1, the partition sum 
is given by 

Zjv = N/[{q)^f . (55) 

Here the factor is a result of the sum of contributions from the A^ ground states and 
the cubic power comes from the product of excitations at the three domain walls. 

Note that the partition sum is linear and not exponential in A^, meaning that the 
free energy is not extensive. This refiects the fact that excitations are localised near the 



Anomalous coarsening and glassy dynamics 



22 



domain walls. In turns this stems from the fact that the energies of most configurations 
are (9(A^^) which is a result of the long-range interaction in the energy function. 

A consequence of this is that in the steady state the system is fully phase separated, 
that is, each of the domains is pure. This was demonstrated in |^ by showing that 

{AA^r) = l-0{r/N). (56) 

for any given r and sufficiently large N. Thus the probability of finding a particle a 
large distance inside a domain of particles of another type is vanishingly small in the 
thermodynamic limit. 

For q close to 1, (g)oo as defined in (^3]), has an essential singularity 



(g)oo = exp [vrVe + 0{1 - g)] | . (57) 

This suggests that extensivity of the free energy could be restored in the double 
limit q 1 and oo with A^lng finite. Physically one can understand this 

scaling variable as the ratio of the domain length (A^/3) to the domain wall width 
(~ jlq^dl/ J g'd/ = 1/1 lng|). The validity of In g as a scaling variable was investigated 
in where a good scaling collapse was obtained for the two point correlation functions. 



3.6. Coarsening 

The analytic results of the previous subsection for A^^ = Nb = Nc give proof of the 
coarsening into three pure domains in that special case. Clearly one expects the same 
behaviour in the general case but to demonstrate it numerically requires prohibitively 
long time scales. In order to study anomalously slow coarsening dynamics numerically 
one can employ an effective or toy model that may be more easily simulated. In this 
subsection I will outline how this can be implemented. 

We consider a system at time t such that the average domain size, (/), is much 
larger than the domain wall width. At these time scales, the domain walls can be taken 
as sharp and we may consider only events which modify the size of domains. This means 
that the dynamics of the system can be approximated by considering only the movement 
of particles between neighbouring domains of the same species. Thus only metastable 
states are considered in the toy model. 

We represent a configuration by a sequence of integers of the form 
ai6iCia2&2C2 • • • cta&aCa where, for example, the zth domain of A particles is of length 
Oj. At each time step a pair of neighbouring domains of the same species of particle, 
say Qi and Oj+i, is chosen randomly. The exchange of particles between domains, takes 
place at a rate dictated by the size of the domains hi and Cj which separate them. Thus 
the lengths of the chosen domains are modified by carrying out one of the following 
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processes: 

1) ai ^ ai-1 



with rate q''^ 

(58) 

with rate q'^' . 



2) Oj — > aj + 1 

Oj+i — s> Oj+i — 1 

If becomes zero, one deletes the domain Oj from the list of domains, and merges bi 
and Cj with and Cj_i, respectively. 

To simulate the toy model efficiently, an algorithm suitable for rare event dynamics 



must be used due to the small rate of events |3^. In an algorithm was employed 
that entails repeating the following steps: 

(i) List all possible events {n} and assign to them rates {r„} according to the rules of 
the model. 

(ii) Choose an event m with probability Vm/R where R = Yln'^^n- 

(iii) Advance time by t — > t + r, where r = 1/R . 

The algorithm would be equivalent to a usual Monte Carlo simulation, where one time 
step is equivalent to one Monte Carlo sweep, if in step 3, r were to be drawn from a 
Poisson distribution R exp[— /2r]. However, a saving in computer time can be had by 
making the approximation t = 1/R. 



In ||2^ the dynamics were simulated for lattice sizes up to 9000. For simplicity we 
consider the case Na = Nb = Nq- An example of typical behaviour of the average 
domain size is shown in figure ^ One can see that after an initial transient growth time 
the data fits very well with a In t behaviour. (Note that the system size is large enough 
that the growth is independent.) Simulations for different q values indicate that, 

(/) = alnt/|lng| (59) 

with a ~ 2.6. The toy model enables one to verify the scaling behaviour (p^ ) and 
estimate the constant a. This would be very difficult to do by simulation of the full 
model (H). 

One can further analyse the toy model by using an Independent Intervals 
Approximation as discussed in Section |1.3|. This was carried out in l2^ . 



3. 7. Generalisation to M > 3 species 

We now generalise the ABC model to M species where M > 3. We may define the most 
general M species model with nearest neighbour particle exchanges that conserve the 
number of each species as follows. Let Xi = 1, 2, . . . , M denote which type of particle 
is present at site i: Xi = m means that site i is occupied by a particle of type m. 
Nearest-neighbour exchanges occur with the following rates: 
g(m, n) 

mn ( ' nm , (60) 
q{n, m) 



Anomalous coarsening and glassy dynamics 



24 



<i> 



100.0 



50.0 



0.0 



10' 



10' 



10' 



Figure 7. Monte Carlo simulation results for the toy model for the average domain 
size, (l), vs. time, t, for N = 9000 and q — 0.8. The data is averaged over 1760 runs. 
There is clear evidence that (l) grows as Int. (Figure taken from l23].) 



and we take q{m, m) = 1. The model conserves Nm, the number of particles of type m, 
for all m. 

According to the choice of the rates (0) the model may or may not phase separate. 
It is not difficult to choose rates so that phase separation does indeed occur. For example 
for M = 4 the choice 

q 



AB 



BC 



CD 



1^ 
1 

q_ 
1 

1 



BA 



CB 



DC 



DA 



AC 



DB 



1 

(1 

1 

_i_ 
' 1 



AD 



CA 



BD 



(61) 



leads to phase separation into pure domains ordered ABCD. Generally, for M > 3 
the structure of the metastable states can become quite complicated |2^. For example, 
domains ordered ACDABCB are also metastable in the model defined by (0). We 
now find the conditions under which the dynamics ( pUD satisfy detailed balance with 
respect to a steady-state weight analogous to ( p] , ^ ) : 

N-l N 

W{{Xi}) = const. X J] J] q{X,,X{) , (62) 



1=1 i=i+l 



where the constant is arbitrary. 
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Consider a particle exchange between sites k and k + 1, where 
and k N. Expanding the product in (|6^), it is easy to verify that 

W{Xi,...,m,n,...,XN) q{n,m) 



n 



(63) 



W{Xi, . . . , n, m, . . . , Xn) q{m, n) ' 
Since this holds for any m, n, and is irrespective of the number of particles of each 
species, the dynamics ( |60D satisfy detailed balance with respect to the weight (|62|) for 
all nearest-neighbour exchanges between sites k and k+1 with k ^ N. If the weights (^) 
are translationally invariant then detailed balance will also hold for exchanges between 
sites 1 and A^. 

Thus, to complete the proof of detailed balance it is sufficient to demand that (^) 
is translationally invariant. To do this we relabel sites i i + 1. The weight then 
becomes 

N-l N 

mm) = n n , (64) 

i=l j=i+l 

where Xq is identical to Xj^. Rewriting this equation by relabelling the indices we 
obtain, 



wax,}) 



N-l N 

n n ^(^^^^ 

i=l j=i+l 



N-l 

n 



N) 



Comparing 



k=l 'i^-^N-'-^k 

with ( p^ ) and noting for example that, 

N M 

\[q{X,,XN) = \[[q{l,X^)f^ , 

j=i 1=1 

one can see that (^2|) is translationally invariant if 



(65) 



(66) 



M r 



n 

1=1 



Q{m,l) 
q{l,m) 



iN, 



(67) 



for every m = 1, . . . , M. Thus, detailed balance holds if ( pTj ) is satisfied. 

In particular for the 3 species model with particles labelled ABC the condition 
(BTI) becomes 



qiA,B)^ 


Nb 




Nc 


\q{B,Ay 


Na 




Nc 


q{B,A)\ 




Uc,A)\ 




[q{A,B)_ 




UiC,B)\ 












\q{c,A^ 


Na 




Nb 














Ub,c)\ 





The ABC model defined in Sec. 3.2 has 



q{B,A) ^ q{C,B) ^ q{A,C) 
q{A,B) q{B,C) q{C,A) 



(68) 
q in which 



case ( |6g| ) reduces to Na = Nb = Nc- 

The ABC model has been generalised to two dimensions |^ where the coarsening 
process generates striped domains perpendicular to the direction of the drive. The Int 
growth of these domains is retained. 
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3. 8. Other related Models 

A model closely related to the ABC model is that first introduced by Lahiri and 



Ramaswamy in the context of sedimenting colloidal crystals. This model comprises 
two interpenetrating sublattices. On each sublattice there reside two species of particles 
with no holes. On the first sublattice (at integer sites) the particles are denoted by + 
and — whereas on the second sub-lattice (at half-integer sites) the particles are denoted 
as 'tilts' \ and /, for reasons that will become apparent below. The dynamics are defined 
as 

ri r2 

+ \- ^ -\+ +/- ^ -/ + 
r2 ri 

(69) 

P2 Pi 

/-\ ^ \-/ / + \ ^ \ + / 

Pi P2 

i.e. the particles make nearest neighbour exchanges on their respective sublattices but 
the rates are influenced by what kind of particle is occupying the intermediate site on 
the other sublattice. In more general rates than (^) were originally considered but 



the above dynamics appear to capture all the generic behaviour 

In the regime ri > r2 and pi > p2 phase separation is observed into ordered domains 
of \— , \+, /+, /— . That is, on each lattice the particles ultimately separate into pure 
domains with one ring rotated by 7r/2 with respect to the other. During the coarsening 
process domain lengths grow logarithmically in time. In |^ the separation into pure 



domains, as first analysed in the ABC model, was termed 'strong phase separation'. 
In the special case where each lattice is half-filled with particles and 

^2 P2 _ .„„v 

— = — = q 70 

n pi 

(where the parameter q has been introduced to allow comparison with the ABC model), 
the steady state weight satisfies detailed balance with respect to a long-range energy 
function very similar in form to (|i2|) 

Wn = Z^\^ (71) 

where 

^ N k 

^= 25ZE^^-V2^- (72) 
i=i j=i 

where r and a are each Ising spins variables for one of the sublattices: cr = 1(— 1) 
corresponds to a + (— ) particle and r = 1(— 1) corresponds to tilt / ( \ ). 

A very useful intuitive picture of this energy function is to regard X]j=i'0-i/2 
the height of an interface (relative to some origin) |]42[| . Thus the dynamics correspond to 
the -|- and — particles making nearest neighbour interchanges on a landscape (implied 
by the tilt variables) that at the same time is evolving in a way coupled to the +, — 
particles. The long-range energy corresponds to the gravitational potential energy of the 
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+ particles on this landscape. A + particle will have its energy minimised at the bottom 
of a valley of the interface. A ground state then corresponds to the configuration that 
allows the + particles to minimise collectively their energy, namely, the tilt variables 
form one deep valley with the + particles residing at the bottom: 

-\-\-\-\ + \ + \ + \ + / + / + / + /-/-/-/ 

The picture of particles moving on an evolving landscape is very appealing in the context 
of glassy dynamics. It is evocative of the glassy regime of a hard-sphere colloid where 
particles move in cages formed by the other particles — when a particle escapes from its 
cage it will cause the cages of other particles to be restructured. 

In the regime ri > r2 and pi < p2 the system is in a disordered phase. In the 
regime ri > r2 and pi = p2 the fluctuations in the landscape are uncoupled to the 
particle dynamics, yet the system still exhibits an interesting coarsening dynamics 



A further model related to the ABC model has been studied by Arndt, Heinzel and 



Rittenberg [^. It was originally couched in terms of +,— particles and holes diffusing 
on a one-dimensional periodic lattice with hop rates 

+ -^^ - + 
1 

+ ^ + (73) 
0-^-0 

In order to make a comparison with the ABC model we identify a + particle with an 
y4, a — particle with a B, and a hole with a C then the dynamics ([73|) become 

q 

AB ^ BA 
1 

CB BC (74) 

AC CA 

Thus this model corresponds to the ABC model with some exchanges forbidden and 
for g < 1 one has the same strong phase separation. However for g > 1 this model 
enters a disordered phase, whereas in the ABC model one has phase separation but 
with the order of domains permuted to ACB. Originally it was thought that there was 
also a 'mixed' phase in the model [^. It now appears that this is a very strong finite 
size effect for g ^ 1 An appealing feature of this model is that the steady state 

can, in principle, be computed exactly by a matrix product approach for all numbers of 



particles ||45 |. 



It should also be noted that a model with cyclic symmetry and non-conserving 
dynamics that exhibits coarsening has been studied |^ . 
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4. Conclusions 

In this paper I have reviewed a variety of simple models that exhibit anomalous 
coarsening in the sense that the dynamics slow down with time and the coarsening 
becomes anomalously slow — slower than the usual power law growth of domain size 
with time. In the models discussed the slowdown in the coarsening is due to dynamical 
constraints rather than any quenched disorder. The dynamical constraints imply that 
energy barriers must be surmounted in order for the system to coarsen further. Thus the 
system is delayed in metastable states for increasingly long times during the coarsening 
process and in this sense the evolution of the system is glassy. 

In Sections 2 and 3 I focussed on two models: the ACIC and the ABC model. In 
the ACIC model the energy barriers encountered during coarsening are logarithmic in 
the domain size whereas for the ABC model the barriers are linear. This leads to the 
domain growth i ~ t'^^"^ (as T ^ 0) for the ACIC and £ ~ Int for the ABC model. 
Despite similarities there are distinctions between the two models. For the ACIC the 
energy function is trivial — free spins in a field — and the dynamics obey detailed balance. 
But it is the existence of forbidden spin-flips that generates the energy barriers in the 
relaxation paths. On the other hand in the ABC model the dynamics are prescribed 
without regard to an energy function and no exchanges (that conserve the particle 
numbers) are forbidden. However in a special case one can identify an effective long- 
range energy function and the anomalous coarsening can be explained in terms of linear 
energy barriers in the coarsening process. Lastly, the ABC model will ultimately coarsen 
into pure domains whose size is related to the system size whereas the ACIC only 
coarsens up to a length 1/e that is independent of system size. For both models there 
exist other related models that exhibit similar behaviour, thus reinforcing their interest. 

The question remains as to how faithfully the anomalous coarsening scenario 
describes a true (experimental) glassy system. Although the models are by no means 
meant to represent any particular system, one recovers the correct phenomenology — for 
example in the ACIC we derived stretched exponential decay of an auto-correlation 
function and KITS law for equilibration. Nevertheless it is widely thought that a 
coarsening system is distinct from a glassy system. 

There are several reasons for this. Firstly, in a glassy system such as hard- 
sphere colloid, the ordered phase (crystalline state) is thought to be irrelevant since 
its nucleation barrier is too large. Thus there are no coarsening domains of an ordered 
crystalline phase. On the other hand at present we do not have any known spatial order 
parameter for a glassy phase so it is unclear how domains of the glassy phase would 
coarsen. Secondly, coarsening implies a definite direction in the dynamics, towards 
the fully coarsened state. Contrary to this, phenomenological trap models of glassy 
dynamics rely on a random exploration of the traps that exist in the phase space. 
That is, when a system manages to escape from one trap it falls randomly into another 
trap. In this way, as the time scale increases, the system will locate deeper traps and stay 
in them for longer. Finally dynamically constrained models by nature rely on specific 
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dynamics rather than the nature of an energy function (if it exists). This contrasts with 
the 'inherent structures' approach to glassy dynamics where the energy landscape is 
the key feature |48, 49]. This point is examined in the present volume by Crisanti and 
Ritort. 

Indeed a criterion has been proposed to distinguish between coarsening dynamics 
and glassy dynamics |pD|] in a microscopic model. One runs a simulation for a certain 
time then makes two copies of the simulation. These copies start from identical initial 
conditions (where the initial simulation was halted) but use different realisations of 
the noise in their dynamics i.e. different sets of random numbers in a Monte-Carlo 
simulation. Then, if the states of the two simulations remain strongly correlated (have 
a finite overlap) as time goes by one has a coarsening system. This is termed type I 
behaviour. On the other hand, if the states of the two simulations become less and less 
correlated then the system is a glassy system this is termed type II behaviour. The idea 
behind this relies on the belief that for a coarsening system there is always a preferred 
direction in the phase space along which all simulations will be swept, whereas in a 
glassy system the traps in the phase space are essentially explored randomly by each 
simulation. 

At present it is not clear how general this categorisation is. For example it is not yet 
clear how systems that coarsen but have many ground states, such as the ABC model 
or the model of E2i, are accounted for. 



Let us also mention a 3-d ferromagnetic Ising model with plaquette interactions ||51||. 
In this model one has the usual doubly degenerate ferromagnetic ground states. In 
addition, flipping any plane of spins does not increase the energy. Thus in total there 
is an exponential number of degenerate ground states. In a quench to low temperature 
type II glassy behaviour is exhibited ||52[| yet it is thought that coarsening occurs whereby 
a characteristic length grows anomalously slowly, possibly as Int. For a related model 
of competing ferromagnet and next nearest neighbour antiferromagnetic interactions 
that has just the two ferromagnetic ground states, logarithmic coarsening is observed 
and explained in terms of energy barriers The glassy coarsening dynamics of the 
plaquette model is not so well understood ||5^. It would certainly be of interest to 
broaden our understanding by studying further examples of anomalous coarsening. 
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